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Abstract 

Recent attempts to recover the graviton propagator from spin foam models in- 
volve the use of a boundary quantum state peaked on a classical geometry. The 
question arises whether beyond the case of a single simplex this suffices for peaking 
the interior geometry in a semiclassical configuration. In this paper we explore this 
issue in the context of quantum Regge calculus with a general triangulation. Via 
a stationary phase approximation, we show that the boundary state succeeds in 
peaking the interior in the appropriate configuration, and that boundary correla- 
tions can be computed order by order in an asymptotic expansion. Further, we 
show that if we replace at each simplex the exponential of the Regge action by its 
cosine - as expected from the semiclassical limit of spin foam models - then the 
contribution from the sign-reversed terms is suppressed in the semiclassical regime 
and the results match those of conventional Regge calculus. 

PACS: 04.60.Pp; 04.60.Nc 



1 Introduction 

Spin foam models jU [2j [31 0] are a discrete, combinatorial and background-independent 
framework for developing a quantum theory of gravity as a sum over histories. A main 
challenge for this approach is to demonstrate that the theory has a well-defined semiclas- 
sical regime in which results agree with conventional perturbative theory [HIE]. Facing 
this challenge involves two different kinds of problems: first, the conceptual problems 
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of defining what is understood by semiclassical limit, and how to connect the different 
languages and conceptual structures of background-independent and perturbative the- 
ories establishing some kind of "dictionary" between them; second, the computational 
problem of finding whether explicit calculations of the same semiclassical quantities on 
both sides of the dictionary really match. 

The last few years have seen great progress on both of these fronts [7|. The starting 
point for this progress was the suggestion by Rovelli [8], [9] of a procedure for comput- 
ing the graviton propagator from Loop Quantum Gravity (LQG) [TJ [10l [11] with the 
dynamics implemented covariantly in terms of a spin foam model. On the calculational 
side, the key ingredients are a boundary semiclassical spin network state peaked on large 
spins and an analytic expression for the large spin asymptotics of the spin foam vertex 
amplitude [12]; on the conceptual side, the framework is the boundary state formalism 
discussed in [1] and in [13], [HI [151 EE E] > which prescribes how to compute observables 
in the boundary of a spacetime region with a path integral over the interior region only. 
If Ox, 2 are local boundary geometry observables (such dihedral angles, 3- 

volumes or lengths [T5], Q3JJ [201 [2D [221 [23] ) acting on a space of spin networks s , then the 
expectation value for their correlation in a boundary geometry q is given in the boundary 
state formalism by 

where ^ q [s] is the boundary quantum state, a spin network functional peaked on the 
classical boundary configuration q = (h, K) (intrinsic and extrinsic geometry), and W[s] 
is the boundary functional, representing a sum over interior spin foams F bounded by s\ 

F,dF=s f e v 

Here Af, A e and A v are the amplitudes that our spin foam model assigns respectively 
to the faces, edges and vertices of the spin foam, which depend on the spins jf and the 
intertwiners i e that label faces and edges for a given F. The graviton propagator can be 
extracted from the connected part of ([I]). 

This expression has been thoroughly analyzed for the Barrett-Crane [21] spin foam 
model P [251 [2E1 [23 [2S1 [21] and has motivated as well the recent introduction of new 
models [3U1 ED E2J E31 El] • However, the analysis has been so far restricted to the case 
where the boundary state has support only on a single graph corresponds to the boundary 
of a single simplex, with the internal spin foam consisting only in the vertex amplitude 
that is dual to this simplex. This is an unsatisfactory state of affairs because then the 
only spin foam variables that are summed over are the boundary spins and intertwiners, 
which are explicitly peaked by the boundary state on their classical values. Therefore 
it could be argued that the nonperturbative nature of the theory plays no role in the 
calculation. A procedure that started with a more general triangulation with many 
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internal vertices, and summed over internal variables in a truly nonperturbative way, 
with the semiclassicality being enforced only at the boundary variables by the boundary 
state, would be a much more nontrivial check of the ability of the nonperturbative theory 
to recover semiclassical results. The question is, then: is semiclassicality at the boundary 
enough to enforce semiclassicality at the interior for a general triangulation? 

A second problem, which could be viewed as a particular aspect of the first one, is 
the "cosine problem": it is known that the large spin limit of the Barrett-Crane vertex 
amplitude is related to the cosine of the Regge action for the corresponding simplex 
[12] , while the conventional Regge path integral involves only the positive exponential of 
the actioiJl]. It was argued in [HI [9] in the context of a single simplex that the negative 
exponential does not contribute to the result because it does not match the phase of the 
boundary state, and its contribution is suppressed. But can something similar occur for 
the case of a general triangulation, where there are many internal simplices contributing 
each a similar cosine factor? 

These are the questions we set out to investigate in this paper, using Regge calculus 
[3U] I3TI ESI EH] as a testing ground. The Regge equivalent of the spin foam boundary 
observable correlation ([!]), for a given arbitrary triangulation, is: 



JUidLtWiLt^Lt) 



(3) 



where now the observables Oi, 2 are functions of boundary edge lengths Lj, ^ q (Li) is a 
boundary state peaking the boundary variables on a classical simplicial geometry q, and 
W(Li) is given by the integral over internal edge lengths of a product over simplices. 



W(L t ) = J X[dl n \[W a (l^L° 



(4) 



W a {l° n ,L1) = ii*{l°n,IA) exp 



h 



(5) 



Here a labels internal simplices, l n are internal edge lengths and (1°, LI) are those internal 
and (if any) boundary edges that belong to simplex a, S a is the Regge action for simplex 
a, and /^(Z^Lf) is the contribution of simplex a to the integration measure \i over 
internal and external edges. 

The goals of this paper are twofold. Firstly, we analyze the semiclassical regime 
of expression ([3]) with a stationary phase asymptotic expansion, and we show that the 
leading contribution to the integral comes from the configuration where all the internal 
variables take the correct classical values (i.e. the particular internal configuration that 
solves the equations of motion for the chosen boundary conditions). The expansion is 



lr This is analogous to what happens in three dimensions for the Ponzano- Regge model |35j 
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first done formally as an h — > limit, but afterwards its regime of physical validity 
is identified with a condition on the boundary state. This is a confirmation that the 
boundary state formalism can enforce semiclassicality for a general triangulation. This 
is a key hypothesis of the graviton propagator calculations that is difficult to investigate 
in the spin foam setting. Regge calculus, as a better-defined and understood model for 
quantum gravity that shares many structural similarities with spin foams, provides an 
ideal arena for testing and checking the validity of this hypothesis. 

Secondly, recall that Regge calculus is actually conjectured to be intimately related to 
the semiclassical regime of spin foam models, via the "cosine problem" mentioned above. 
Accordingly, we investigate the consequences of replacing (jHJ) by a similar expression 
that is conjectured to result from the large spin limit of the spin foam vertex amplitude 
corresponding to each simplex: 

W a {r n) L°) = LI) [P*(l°, LI) cos (S CT (C LI) + |) + D a (l°, Lf )] , (6) 

where the term D a is dominant but slowly varying; this ansatz is motivated by the re- 
sults on the asymptotics of the Barrett-Crane model [121 EH El H2]- We show that, 
remarkably, the results of this replacement when the expression is introduced in (j3J) are 
the same as those obtained from the original Regge expression: the differences between 
(12]) and (jHJ) represent deeply quantum configurations that drop out from the semiclas- 
sical regime of boundary observables. This is an important step towards showing that 
the correct semiclassical limit can be indeed recovered from spin foam models, making 
contact with perturbative quantum gravity a less unmanageable goal. 

The paper is organized as follows: Section [2] is dedicated to the purely Regge part of 
the calculation. Subsection 12 . II introduces the classical Regge action, the path integral for 
the quantum theory, and the semiclassical boundary state we will consider. Subsection 
12.21 contains the explicit evaluation of boundary observables in Regge calculus, showing 
with a stationary phase analysis that the Regge path integral peaks on a single semiclas- 
sical configuration. Subsection 12.31 discusses the particular case of flat space and how to 
deal with the translational invariance that arises. Then Section [3] repeats the main cal- 
culation for the modified path integral resulting from showing that the semiclassical 
results are unchanged and therefore that, under certain assumptions, the semiclassical 
regime of spin foams matches that of Regge calculus. Section 0] is a discussion of the 
results and suggestions for further work. 

Throughout this paper we work in four- dimensional space with Euclidean signature. 
We use units with c = 1 but keep h and k = 8ttG explicitly in order to keep track of 
semiclassical expansions. The Planck length is then lp = s/Hk. 
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2 Semiclassical regime of Regge calculus 



2.1 Regge calculus and boundary state 

The Regge action for Riemannian four-dimensional simplicial gravity [36l H3j HU H5] can 
be written as 

5 = -VU. (7) 

K / — ' 

a 

Here a ranges over the faces of the triangulation, with A a representing the area and 4> a 
the deficit angle at each fac^E The deficit angles are defined by 

a = 27r-^0;, (8) 

where 9° is the internal dihedral angle at face a of simplex a, and the sum is over all 
simplices meeting at face a. If the manifold has a boundary, the action takes the same 
form with the deficit angles for boundary faces given instead by 



7T 



£*.*■ (9) 



The action ([7]) is to be considered a function of the edge lengths l iy with i ranging 
over all edges. The Regge equations of motion obtained from the variation of the action 
are: 

= < 10 > 

a 1 

The quantum theory is normally defined, for a given triangulation, by the path inte- 
gral over all edge lengths in the triangulatiorjfl [37J EH [39] : 

Z = J Hdkvik) exp (jrSik)) , (11) 

where /x(Zj) is a suitable measure function [17J HHJ 0H [50] and some fixed asymptotic 
boundary condition is understood. Following the framework outlined in the Introduction, 
we use instead the boundary formalism in which the lengths integrated over are only those 
in the boundary (Lj) and the interior (/„) of a closed, compact region, with a boundary 
quantum state ^f q (Li) included in the path integral: 

Z q = j \[ei n n{l n ,Li) exp (^S{l n , Li^j %(Li) . (12) 



2 Most of our results would be unchanged if we included in the action a cosmological term summing 
over the volume of simplices, — ^ Vo-- We do not include explicitly this term for the sake of brevity, 
but we will point out any point where it would make a difference. 

3 "Problem 57" in \M- 
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Expectation values of boundary observables O acting on the space of boundary states 
are denned by 

(d) q = y Jll dLi H^nKln,^) exp^S^L*)) 6%(Li) . (13) 

The choice of boundary state can be taken as equivalent to the choice of asymptotic 
boundary conditions for the full path integral; in fact, as we discuss in Section HI the 
boundary state can be defined as the result of integrating out the external edge lengths 
in fjTTTT pl. We refer to [51] for a discussion of this procedure in the case of linearized Regge 
calculus on flat space. Aside from the discussion in Section HJ in this paper we will not 
be concerned with the first-principles derivation of the boundary state. We take instead 
the "phenomenological" route that assumes the boundary state to be a Gaussian in 
the fluctuations (plus eventual higher-order corrections) around some specified classical 
boundary configuration {L Q m }. This follows the approach used in [SJ [H [251 [2Z1 [2H1 [2H] for 
spin foam models and in [27] for Regge calculus with a single simplex. For a discussion of 
the role of the boundary state in the three-dimensional case with a single Ponzano- Regge 
vertex see [521 [531 M, ES] . 

More concretely, the recipe we use for the boundary state is as follows. Given a 
triangulatiorjfl and the 3d subset of it that is the boundary, we want the boundary state 
to be a function of the boundary lengths {Li} that selects a particular semiclassical 
regime for the quantum theory: one that is peaked on a particular classical solution. 
For this we first specify a set of values for the boundary lengths, Lj = L\. When used 
as Dirichlet boundary conditions in the equations of motion fllOp . there may be for the 
internal lengths either (a) no solution, (b) a unique solution, (c) a discrete set of solutions, 
or (d) a continuum of solutions. We are not interested in the case (a), because our goal 
is to write a state that peaks us on a classical solution; let us therefore assume that we 
choose {L®} so there is at least one solution. 

Focusing on case (b) first, call the unique internal solution l n = Take the 

derivative of the internal action with respect to the boundary edge lengths, and evaluate 
it at the classical solution: 

dS 



(14) 



I -1° 



This is the boundary extrinsic curvature of our chosen configuration. An alternative 
definition of is the derivative of the Hamilton function with respect to the boundary 
length Li. The Hamilton function S H is the internal action evaluated at the classical 



4 A key property making this viable is that the action S and the measure fj, are both local, splitting 
respectively into a sum and a product over simplices of single-simplex actions and measures. Therefore 
talking about the "internal" action and measure of a region is well-defined. 

5 Throughout this paper by triangulation wc mean only the topological triangulation, also known as 
the skeleton. 
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solution, as a function of boundary variables, and the equivalence of both definitions 
comes from the expression 

dS^ = d^ + d^dln 

dLi dh i dl n dLi ' 

where the last derivative is evaluated through the dependence of internal variables on 
external ones once equations of motion are solved. When evaluated at a classical solution 
the second term vanishes and the equality is established. 

The Gaussian boundary state that we use is then given by: 



V q (Li) = C exp 



2hn 



a l3 (Li - L°) (Lj - L 



0\ 



exp 



(16) 



where the matrix otij is symmetric and has positive-definite real part, C is a normalization 
constant, and summation over i,j indexes running over boundary edges is understood. 
The reason the phase of the Gaussian involves K® is that Ki is the canonically conjugate 
variable to Lf, the state *f? q is then peaked in both "position" and "momentum", as a 
good semiclassical state should. See [HI IH] for a more extended discussion, and the next 
subsection of this paper for a discussion of the consequences of using the wrong phase. 

If we are in case (c), then we need only a slight modification to the above procedure. 
The boundary lengths Lj = L° are not enough to specify a solution. To get a state 
peaked on a particular solution of the many compatible with the Dirichlet condition, we 
choose the solution we are interested in, calculate its extrinsic boundary curvature Kf 
as above, and construct the Gaussian boundary state (fT6|) in the same way. Then this 
boundary state selects the semiclassical regime peaked around the particular solution 
chosen. 

Case (d), finally, is what happens when one solution compatible with the given bound- 
ary lengths includes a patch of flat space: a region where one or more internal vertices 
of the triangulation can be translated keeping the physical geometry invariant, thus 
producing for some of the edge lengths a continuum of values representing the same 
simplicial geometry. In this case, the action evaluated in all these solutions will be the 
same - the Hamilton function of the configuration - and therefore there is a well-defined 
extrinsic boundary curvature K® given by its derivative with respect to the boundary 
length L°. The corresponding Gaussian boundary state peaks us on the single interior 
geometry compatible with the given boundary conditions, despite not fixing a unique set 
of values for the internal edge lengths. We will see in Section 12.31 that this is enough for 
a well-defined evaluation of observables. 

To conclude this section, note that if the coefficients are scale-independent (i.e. 
they are invariant under a global rescaling of lengths) then the Gaussian fluctuations 
for both conjugate variables vanish asymptotically for L ^> \^hn, where L is a typical 
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length scale of the chosen classical solution: 



ALi AKi VHk r— , , 

-jjr~l^r~— >°> ^o»v^. (17) 

Thus our state is sharply peaked on a classical configuration, with small relative disper- 
sion, when the typical length scale of this configuration is much larger than the Planck 
length. The ratio of the Planck scale to the boundary scale will serve as a physical 
expansion parameter in our calculation, and we will evaluate only the lowest nontrivial 
order of observables in this expansion. This justifies us in using a Gaussian boundary 
state; higher-order corrections can be taken into account order by order in the asymptotic 
expansion using an improved boundary state, as discussed in [27] . 



2.2 Stationary phase evaluation of observables 



In this section we show explicitly how the computation of boundary geometry observables 
proceeds within our framework. To the order of accuracy allowed by restricting the 
boundary state to the Gaussian form (fl6j) . a boundary observable is given by: 



(6) q = — I Y[^ndLifi(l n ,Li) exp ( ^S(l n ,Li 



O 



x exp 



2Hk 



aiij (Li - L°) (Lj - L 



exp 



Here O may be an intrinsic boundary geometry observable, i.e. a function of bound- 
ary lengths Li that acts simply by multiplication, or an extrinsic boundary geometry 
observable, i.e. a function of the boundary extrinsic curvature Ki that acts on ^> q as a 
derivative with respect to Lj. 

Expression ffl8l) is a well-defined integral, which could in principle be computed nu- 
merically (after a prescription for the measure is chosen). We are interested in evalu- 
ating it analytically in a semiclassical approximation. The idea is that the boundary 
state enforces the integral to be peaked in a particular classical internal configuration. A 
"picture" for how this happens can be seen by doing the integral layer by layer (cf. |56j): 
starting by integrating out the boundary lengths that are explicitly peaked on classi- 
cal values, and then moving towards the interior, with each "layer" of edges getting 
peaked by the immediately previous one in the configuration that solves the inwardly 
propagating equations of motion. 

Carrying out this procedure explicitly may be feasible for particular examples, but 
it is complicate and makes difficult the extraction of information about boundary cor- 
relations. A much simpler technique that works for general configurations is to use a 
global stationary phase approximation. This can be done formally as an h — * limit, 
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or in other words an asymptotic expansion of the integral for large values of hr 1 . We 
will show that this peaks the whole integral in the global configuration that solves the 
equations of motion for the given boundary. An a posteriori examination of the result 
will exhibit that this formal asymptotic expansion provides in fact a good approximation 
if a physical condition is satisfied: the ratio of the typical scale selected by ^ q to the 
Planck scale must be large. The small h stationary phase approximation is a way to 
shortcut to an asymptotic expansion in this physical parameter. 

Examining ffl8|) we see that if the argument of the exponentials was a purely imaginary 
quantity we could use a straightforward stationary phase approximation, while purely 
real (and negative) argument would require a Laplace approximation. But the argument 
in ( [181) is complex, so we need a more careful statement of our approximation. Let 
us therefore make a brief mathematical detour to state the precise result that we will 
apply and also the formulas for evaluating successive orders in the asymptotic expansion, 
beyond the lowest-order expression that can be found in textbooks. 

We are interested in estimating an integral of the form 

„ d 

F(X)= / Udx' f(x) e lX ^ x ) (19) 
Jd i=i 

for large parameter A. In (1191) . D is an open region in l d , / and Q are two smooth 
functions / : D — > K and Q : D — > C. We assume that / vanishes on the boundary of 
D. Let Hij(x) be the Hessian of Q(x): 



We assume that the imaginary part of the Hessian is positive semi-definite, lm[Hij(x)] > 
Vx G D. As a result of the eventual presence of zero eigenvalues of lm[Hij(x)], we have 
that the set of minima of Im[Q(x)] is a subset of D, which we call Dq. Let us assume 
for the moment that Re[Q(x)] has at most one stationary point. There are then three 
possible scenarios: the real part of Q(x) has (i) no stationary point in D, or (ii) it has a 
stationary point x in D, but it does not fall in the subset D , or (iii) there is a stationary 
point x falling in Dq. 

We can state the following extension of the stationary phase approximation JSTJ [55] : 
in the cases (i) and (ii), for large parameter A the function (I19p decreases faster than any 
power of A" 1 , F(X) = o(X~ N ) ViV > 1. On the other hand, in case (iii), the function F(X) 
is of order X~ d ^ 2 and (assuming Det[Hij(x)} ^ 0) for large A the asymptotic expansion 
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of the integral i^| 

/arN* JXQM (f anA - n + o(A - K) V (21) 

1 ' \XJ ^\AeiH,,(x)\ ^ " 'I 1 

In (1211) . the index ind [iJjj] is defined in terms of the eigenvalues Hj{H) as ind[i?y] = 
2 arg|/ij with — | < arg[/ij(if)] < +|. The asymptotic expansion coefficients a n 
can be expressed in terms of derivatives of the functions f(x) and Q(x) evaluated at the 
stationary point x. 

We define R(x) to be the difference between Q(x) and its second-order Taylor poly- 
nomial at x, R(x) = Q(x) — Q(x) — ^Hij(x)(x — x) l (x — Such function has a Taylor 
expansion around x which starts with a term cubic in (x — x) 1 . The coefficient a n is 
given by^ 



(23) 



J=0 



In particular, the zeroth-order and the first order coefficients are given by 

a = f(x) , (24) 

ai= i\fS\x){H~\x)f+ 

+ i \f{x) R^l(x) R%> nl ( x ) {H~\x)r (H-\x)y n (H-\x)f . 

The coefficients a n admit a Feynman-diagrammatic representation in terms of propaga- 
tors (iJ~ 1 (x)) u , vertices Rf^ (x) with b > 3 and insertions corresponding to derivatives 



6 If there are several isolated stationary points with non-degenerate Hessian within Dq, we obtain a 
sum of contributions of the form (|2ip from each of them; if there is in addition one or more stationary 
point outside Do, its contribution is o(\~ N ) ViV > 1 and hence suppressed in comparison with the 
other ones. The only remaining case of interest, which is a continuum of stationary points within Do, 
will be considered in the next section. 

7 The formal expression f(x + -^j) and the analogous one for R(x) stand for a Taylor expansion of 
the functions around the point x with the fluctuation substituted by a derivative 

/(i+ ^ ) = /w+ |J. / <« (i( , ) _|_.._|_. m 

Notice that for fixed n in ([23]) only a finite number of derivatives appear. Here fj^\ (x) denotes the 
&-th partial derivative of / with respect to x% x ■ ■ ■ Xi b , evaluated at x. 
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of / at x. Since only diagrams with n loops contribute to a n , this kind of approximation 
is called "loop expansion" in the context of quantum field theory [59] . 

Now we go back to our physical problem to see how these mathematical formulas 
apply. As we said, fir 1 plays the role of large parameter A for the stationary phase 
evaluation of (IT8"|) . Our domain D consists in a subregion of M M+7V , where M and N are 
the number of boundary and internal edges in our triangulation; let us call points in this 
space x = (Li, ■ ■ -L M ,li, ■ ■ - In)- The subregion is specified by the measure fj,(L m ,l n ), 
which vanishes for negative values of the lengths, and also for values not satisfying a set 
of triangular inequalities that make the configuration possible^! The function we need 
to extremize is 

Q(l n , Li) = S(l n , Li) + ay (Li - L°) (L, - L°) - K° (L t - L°) , (25) 
and its Hessian is given in block matrix form by 

H-( A ;Tb) (26) 

where the block separation corresponds to boundary and interior variables, in that order, 
and 

r » n 9 S r -, OCoH r-« -^n 9 S r _„ -. 9 S , „ _v 

W'^W M « = V- l N l™=alA- l B l- = S^' < 27) 

Our assumption that Im[iJ] is positive semi-definite is clearly satisfied thanks to the 
properties we have stated for a. It is clear that the subset Dq of the domain where Im[Q] 
finds its minimum will be given by points of the form 

x = (l n ,L° i ); (28) 

that is, the condition that the imaginary part of Q be minimized imposes the boundary 
lengths at the extremum to be equal to the classical values specified in the boundary 
state, and it imposes no further condition on the interior lengths. 

Now we look at the real part of Q(x); taking the above condition coming from the 
imaginary part as given, the stationary points will be the solutions to the equations: 

= , (29a) 

^•S* (l n , L t 



-Kf= . (29b) 

L t =L° 



s We assume that fj, vanishes smoothly at dD so there are no boundary contributions to the integral. 



11 



The first of these amounts to the equations of motion with boundary conditions L\. 
Assuming we are in case (b) defined in the previous subsection, there is a single solution 
for the internal lengths, l n = l^(Lf). The second equation can then be simultaneously 
satisfied only if the phase factor of the boundary state, Kf, matches the extrinsic bound- 
ary curvature at the classical solution. A similar thing happens if we are in case (c), 
the difference being that then the second equation is not just a consistency condition 
but actively selects the desired classical internal solution out of the several compatible 
with the boundary lengths. In both cases, if the phase of the boundary state is an arbi- 
trary quantity unrelated to the solution(s) implied by the boundary lengths, we fall into 
case (i) or (ii) of the theorem and the semiclassical limit is suppressed. This confirms 
that only a state properly peaked in intrinsic and extrinsic classical geometry can select 
correctly the semiclassical regime of the theory. We stress that our formalism does not 
merely peak the path integral in any configuration that solves the classical equations of 
motion, but only in a specific classical solution, namely the one that solves the equations 
of motion for the boundary conditions encoded in the boundary state. 

We still need to check that H(x) is invertible. As far as the Hessian of the Regge 
action is concerned, it is known [60l ED] that it has in general a zero mode corresponding 
to global conformal rescalings, and also translational zero modes in the particular case 
of patches of flat space. The conformal zero mode disappears from H when we add 
the iat part coming from the boundary state; physically, this corresponds to the fact 
that the boundary state fixes the geometry of the boundary and in particular its sea 
The translational degeneracy persists, as can be easily checked; this corresponds to 
the fact mentioned (d) in the previous subsection, that if boundary conditions 

corresponding to an internal solution that is totally or partially flat space are fixed, then 
there is a continuum of internal classical solutions, which correspond to translations 
of interior vertices of the triangulation. For the moment we assume that the specific 
solution selected by our boundary state does not contain flat regions, and then there 
will be a single extremum with nondegenerate H. The case of flat space is left for the 
following subsection. 

In summary, the point in configuration space from which the integral gets its dominant 
contributions as specified by the stationary phase approximation is (L°, i°(L°)): the 
boundary lengths at their chosen classical values, the internal ones at the classical values 

9 The conformal zero mode reflects the fact that, given one classical solution, a family of classical 
solutions can be found (in absence of matter and of cosmological constant) simply rescaling all the lengths 
by a given parameter. However, these are new solutions: they represent simplicial geometries with 
observables physically larger or smaller than the original ones. In the quantum theory, the semiclassical 
regime is identified by a semiclassical boundary state, which is peaked on a specific intrinsic boundary 
geometry with a specific scale; this selects one among the classical solutions related by a rescaling, and 
therefore the conformal mode disappears from the Hessian when a is included. A different perspective is 
taken in [51] , where the conformal zero mode is considered a gauge degree of freedom and a gauge- fixing 
procedure is advocated. 
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determined from the boundary conditions via the equations of motion. 

Let us consider next the evaluation of observables. We focus on intrinsic boundary 
geometry observables depending only on boundary lengths, but a parallel discussion 
could be done for extrinsic boundary geometry. Applying our collection of formulas, the 
expectation value of an observable 0(Lj) is given by 

(0(Li)) q = 0{Lfj + order H corrections . (30) 

This comes from keeping only the zeroth order in the asymptotic expansion both for the 
numerator and the denominator Z q ; notice that dependence on the measure /i drops out. 
Note also that, while the "kinematical" expectation value of an observable (the result 
of integrating it against the modulus squared of the Gaussian boundary state over the 
boundary variables Lj) is exactly 0(L®), the "dynamical" expectation value calculated 
with ( lT8i) acquires corrections of order h. 

We are specifically interested in correlations of boundary observables, which are of 
order h. The simplest example is the connected two-point function for lengths, which 
is related to the simplicial graviton propagator studied in [601 [6TJ [62]. Applying our 
asymptotic expansions to the first order gives: 

( U Lj ) q — {Li) q { Lj ) q = h + 0(h 2 ) (31) 

= h [(A + ia — N* B^N)- 1 ]^- + 0{h 2 ) . (32) 

Recall that H is a (M + A)-dimensional matrix involving second derivatives of the action 
with respect to boundary variables Lj and internal variables l m . The correlation function 
involves taking its inverse and selecting from it the M x M block corresponding to the 
boundary indices. According to the well-known formulas for inverting block matrices of 
the form (|2"6"I) . we obtain the second line of (I3"2"l) . 

Expression (1321) has actually a very simple and beautiful meaning. Consider the 
Hamilton function S H (Li) = S(l m (Li), Li): the action evaluated at the classical solution 
for given boundary lengths {Li}, considered as a function only of these boundary lengths 
via the equations of motion. We have that the second derivatives of the Hamilton function 
can be written in terms of derivatives of the action evaluated at the classical solution in 
the following way 

d 2 S H Q _ d 2 S / d 2 S \-i d 2 S d 2 S 
dLidLj 1 dLidLj \dl m dl n ) dLidl m dLjdl. 

This identity can be easily derived using the implicit function theorem, considering the 
first derivatives of S(l m , Li) being set to zero as describing a constraint surface that makes 
the internal variables U an implicit function of the boundary variables Lj. Comparing 
(|33|) with (I32p . we see that the two-point correlation function to lowest order is given by 



(33) 
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the inverse of an M x M matrix which is the sum of two terms: the Hessian matrix of 
the Hamilton function and the correlation matrix of the boundary state 



(L i L j ) q -(L i ) q (L j ) q = h(-—(L°) + -a ij ) . (34) 

Note that since the action includes a prefactor k~ 1 both matrices scale in the same way: 
the two-point function is given by Hk = lp 2 , multiplied by a dimensionless matrix. For 
fixed boundary state, the "scaling with distance" of the two-point function (the fact 
that correlations decrease with distance in a way compatible with the 1/r 2 behaviour of 
continuum gravity) is hidden in the way (1341) changes when considering different pairs 
of edges: edges that are further away correlate more weakly than nearby edges. 
Recovering the exact scaling known from conventional perturbative theory on a regular 
lattice [601 EE] can be imposed as a condition on the boundary state matrix a. This is 
discussed more extensively in Section HI 

From the correlation function for lengths (1341) . all two-point correlation functions of 
local intrinsic boundary geometry observables can be obtained accurately to the lowest 
order. In particular the £1X6 cL~ 3X6 cL correlations, which were computed for the case of a 
single regular simplex in [27], can be easily obtained for a general configuration. As 
is generally the case in semiclassical perturbation theory, the correlation functions are 
to the lowest order independent of the nonperturbative measure /i. In a Feynman- 
diagrammatic language, at tree level the measure does not contribute, but the one-loop 
0(h 2 ) correction does depend on the measure. Finding the general form of the second- 
order result would also involve using an improved boundary state of the kind discussed 
in [27]. 

To conclude this section, let us discuss which is the physical regime in which the first 
terms of our formal asymptotic expansion provide a good approximation for the observ- 
ables. For this we compare the computed expectation value of a boundary observable to 
its classical value, and require that 

Computing this with the stationary phase asymptotic expansion for a general observable 
leads to a condition on the classical solution we have chosen to peak upon. Consider a 
classical solution which oscillates over a length scale A. From the general form of the 
expansion coefficients a n it can be seen that 

f ~ Hk/X 2 . (36) 

This is analogous to the conditions of validity of the WKB approximation in quantum 
mechanics [63]. Recall in turn that the classical solution is uniquely determined (up to 
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the flat space degeneracy discussed in the next section]^ by boundary lengths L\. In 
absence of matter and cosmological constant A (which introduce intrinsic length scales 
in the classical theory), the Regge action is a homogeneous function of the edge lengths. 
Under these conditions, a global rescaling of the edge lengths by a factor L leads to 
an identical rescaling in all the internal lengths at the classical solution determined by 
the new boundary lengths. Hence our parameter X will also scale with L . We conclude 
that we can take Hk/Lq as the physical parameter of our expansion, with L being a 
typical scale of our boundary lengths L\. Thus the semiclassical regime is characterized 
by having boundary lengths very large with respect to Planck length: 

l -f~ 2 « 1 • (37) 

Note that this ratio is precisely the parameter that quantifies the kinematical relative 
dispersion in the boundary Gaussian according to (fTTI) . If the action includes a cos- 
mological term or matter sources, then (1371) is not generally true, while f)36p keeps on 
holding; in this case, careful study of the classical solution is needed to identify a good 
expansion parameter, which is no longer a single global scale of the boundary state. 

2.3 Flat space 

Let us now go back to (fTBj) and consider the case we labelled (d) in Section 12.11 the 
classical internal solution implied by our boundary data includes patches of flat space. It 
is easy to construct such solutions; for example, one can start from any solution defined 
on a given triangulation, and then modify the triangulation adding an extra vertex 
inside a simplex, specifying the lengths of the new five edges to be such that the ten new 
triangles that include these edges all have zero deficit angle (a "one-five" Pachner move). 
See [61] for a recent discussion of such configurations. When the classical solution that 
our boundary data selects includes flat regions in this way, the edge lengths are redundant 
variables, in the sense that there is a continuum of values for them - corresponding to 
translations of some internal vertices of the triangulation - that correspond to the same 
physical geometry. As a result, all solutions for the edge variables within this continuum 
of values have the same value of the action, which is the extremum for the given boundary 
data. 

If the number of internal vertices that can be translated leaving the Hamilton function 
unchanged is q, then the N internal length variables l n can be divided in one set of 4q 
variables parametrizing the position of the points in the flat four- dimensional interior 
regions, and a set of N — 4q variables that are in a sense "orthogonal" to this subspace. 

10 In case (c), where more than one classical solution correspond to given boundary lengths, the 
argument that follows keeps on holding since K®, which determines the solution in conjunction with L° Ll 
scales with Lq as well. 
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For given values of the first Aq variables, there is a unique set of values for the remaining 
N — Aq variables that extremizes the action. Thus it is possible to do first the integrals 
in the M boundary variables and in N — Aq of the internal variables using a well-defined 
stationary phase approximation of the form described in last section: the boundary 
variables get peaked in the values L°, and the N — Aq internal variables get peaked in 
the values that extremize the action for fixed value of the remaining Aq variables^. This 
leaves only the remaining 4g-dimensional integral over all the possible flat configurations, 
which share the same physical geometry. 

Using the notation x n for the N — Aq internal variables integrated over, calling the 
other Aq internal variables y n , and calling x1(y n ) the values of the x n that extremize Q 
for given values of y n (and L°, though we leave this dependence implicit), the result we 
obtain after the first set of integrations is done by stationary phase is: 

(° > = T / ll d ^" ri n „ Tn n , : — rjTTa a A L t x n{Vn ),y n ) • (38) 
Z q J t\ [det H{Ll x» n (y n ),y n )] 1/2 Y 

where Ti stands for the Hessian matrix of second derivatives of Q with respect to the 
variables {Li,x n }. As usual, the coefficients a,j are defined by (|23j) with the role of 
Q being played by (1251) and the role of / by n O. Z q is identical to the numerator but 
omitting O. Note that ( 1551) allows well-defined computation of the boundary observables 
to any desired order; the boundary state peaking on a specific geometry is enough and 
there is no need to "gauge fix" in addition the position of the internal vertices in flat 
regions. 

To the zeroth order, we have as before that in the numerator ao equals the measure 
/i(L°, x^(y n ),y n ) times the insertion C(L°). The second factor does not depend on the y 
variables, and upon taking it out from the integral we are left with a "volume" integral 
over the flat regions that cancels out with Z q . Hence the result of the previous section 

(0(Li)) q = C(L°) + order h corrections . (39) 

is still true when flat regions are taken into account. 

What about the two-point correlation function? Taking (1551) up to 0(h) for this 
observable leads to 

(Li Lj) q — {Li) q {Lj) q — (40) 



h 



I UtLi d Vn [det H(V?, xl(y n ),y n )} 1 2 ^x(L°, x° (y n ), y n ) ?^/(L°, x° n (y n ),y r , 
IUti d Vn ldetH(Llx°(y n ),y n )}- 1/2 ti(Lla*(y n ),y n ) 



The argument used in the previous subsection to relate the Hessian of Q to the Hessian 
of the Hamilton function still applies to 7i, even though it does not include derivatives 
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This is akin to the so-called method of collective coordinates in statistical field theory. 
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with respect to the y n variables, because these are precisely the variables the action is 
independent of. Hence in ( HU1) we can replace TC^ 1 by f ^§^- (L®) +iaij/i^j , which 

is transparently independent of the y n variables. We can therefore take it out from the 
integral in the numerator and then cancel out the remaining "volume" integrals, leaving 
us with 

( U Lj ) q — (Li) q ( Lj ) q = h (^^ra + ~* ■ (41) 

which is the same expression we had in the previous subsection. The conclusion is that 
the first-order correlation function can be found directly from the Hamilton function of 
the configuration, when there are flat regions just as when there are not. 

An interesting implication of this is that the leading order of the correlation func- 
tion in a flat configuration is triangulation-independent. All our calculations so far have 
assumed a given fixed triangulation. But consider now that we choose a different tri- 
angulation of the interior region, one for which the internal solution selected by the 
fixed boundary state is (at the level of physical geometry) the same one as before; the 
differences between them amount only to a series of Pachner moves performed strictly 
within flat regions of the original solution. Since this transformation leaves the Hamilton 
function unchanged^], it also leaves unchanged the first-order quantum correlations of 
boundary geometrical observables. We do not know if this result is preserved at higher 
orders in the asymptotic expansion; in principle, there seems to be no reason to expect 
it. 

These results are of great interest for the ongoing discussion on the meaning and scope 
of triangulation independence in spin foam models. They suggest, for example, that spin 
foam observables computed with the boundary formalism for a configuration with five 
spin foam vertices (which amount to five simplices in the original triangulation) might 
equal those for a single vertex, at least to the first order in a semiclassical expansion 
peaked on a flat configuration. 

3 Semiclassical regime of spin foam models 

In this section we undertake the extension of our results for Regge calculus to the context 
of spin foams. A full discussion of the semiclassical limit of spin foam models would 
require starting with a model and attempting to compute correlations of observables 
on a semiclassical state, as outlined in equations ([I]) and ([2]). We should then show 
that - by some reason analogous to the stationary phase analysis of Regge calculus - 
the contribution coming from the desired internal semiclassical geometry dominates the 

12 This is because the Hamilton function depends only on the classical simplicial geometry and not on 
the triangulation and the assignment of edge lengths separately. 
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state sum, and the results can therefore be related to those of conventional perturbative 
theory. 

Fulfilling this outline in its entirety is beyond the scope of this paper. We focus instead 
on examining one particular problem that might hinder its completion - the "cosine 
problem" mentioned in the Introduction. The conjectured semiclassical limit of the spin 
foam vertex amplitude is not merely the exponential of the Regge action that appears 
in the Regge path integral, as would seem to be needed to get a perfect correspondence 
with simplicial gravity, but instead features the cosine of the Regge action plus a term 
with no clear geometrical interpretation [T2], |40j UH |42] . It was shown in [H [9] that this 
does not thwart the correct semiclassical limit for the case of a single simplex, but doubt 
remained as to whether this success would survive when more general triangulations are 
considered. In this section we address this question within the general framework of 
this paper, replacing at each simplex in the Regge path integral the exponential of the 
action by its spinfoam-inspired equivalent, and showin g th at the contribution from the 
non-Regge terms is suppressed in the semiclassical limit] 13 !. 

Let us first state some expectations for the semiclassical behaviour of the vertex 
amplitude of a spin foam model (for general reviews of the spin foam framework, see 
[21 OH]). A spin foam vertex amplitude A v (jf,i e ) for a four- dimensional model depends 
on the spins jf and the intertwiners i e that colour the ten faces and five edges surrounding 
the vertex; these faces and edges are dual to the triangles and tetrahedra surrounding 
the simplex dual to v in a simplicial triangulation of the manifold. When the spins jf are 
sent to infinity and the intertwiners i e are chosen to be a specific function of the spins, 
the vertex amplitude is conjectured to have the following asymptotic behaviour: 

+ D V , j f -> oo . (42) 

In this expression, P v scales as j^ m when all spins are rescaled by a factor j . The 
coefficients Of do not scale with j Q and each has the appropriate value for the internal 
dihedral angle at the triangle dual to face / in the simplex dual to v. Newton's constant 
enters in the argument of the cosine through the identification of the spins jf with the 
areas of the triangles by the LQG formula A a = Hnjjf for the large spin asymptotics 
of the spectrum of the area [HJ [19]. Note that the argument of the cosine is therefore 
the Regge action for the simplex dual to v (divided by fij^j. The term D v , finally, does 

13 This can be seen as a continuation and generalization of the work done in [SS] for the Ponzano- Regge 
model in 3d. See also [551 EH GUI EH] for a different perspective on the "cosine problem" that traces it 
to the need to implement causality in the model. 

14 Recall that, in order to obtain a deficit angle, we have to subtract to 2tt the sum of the dihedral 
angles at a triangle. As a result, to recover the Regge action from a sum over simplices, an extra term 
2irA a for each triangle is needed. This term is expected to arise from the face amplitude of the spin 
foam model. 



A v (jf, i e ) ~ P^COS 
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not have any "good" interpretation as a classical geometry. It scales as with d < m, 
and is therefore dominant over the cosine term, but it contrasts with it in being slowly 
varying (not oscillating). The presence of a cosine, of the 7r/4 and of the dominant non- 
oscillating D v term are standard features of asymptotic expansions. The form (l4~2l of 
the asymptotic vertex amplitude is known to be true for the Barrett-Crane (BC) model, 
with m = 9/2 and d = 2 [T^l E2 SD H2]- This is the main motivation for us considering 
it, under the conjecture (made reasonable by recent developments discussed in the next 
section) that the more recent spin foam models [30], ED E2] [33], [31] lead to identical or 
similar behaviour. The only addition we have made to the BC result is the Immirzi 
parameter 7 to scale the areas correctly in accordance with LQG; since the new models, 
unlike BC, include 7 and reproduce the LQG area spectrum [31], this seems a reasonable 
assumption. 

We want to focus our investigation on the consequences of having such an expression 
at each simplex of the triangulation. We therefore examine the following definition of 
observables: 



PM; r L-)eos(±s a (i:,i4) + j) 4 dm;,.l:) 



O exp 



2hn 



a>ij (Li - L?) (Lj - L°) 



exp 



(43) 



This is identical to the Regge calculus expression ffl8l) . with the exponential of the 
action replaced, at each simplex a, by the conjectured limit of spin foam models. The 
dependence of S a , P a and D a has been translated from the spins to the lengths via the 
areas; (/£, Lf) are the internal and (if any) boundary lengths belonging to simplex a, 
and we write \x for the measure because it may be different from the pure Regge measure 
\i. There are several hypothesis and assumptions needed for using (H3|) as the correct 
semiclassical regime of spin foam observables; the plausibility of them will be discussed 
and defended in the next section, but for the moment (]4"3l can be taken as an ansatz 
for focusing on the effects of the U P cos S + D v structure and studying whether they 
hamper the correct semiclassical limit. 

The product over simplices can be expanded using Euler's formula leading to 



(6) = Yj H dl ndLmKln,Li) ^ 



HP:(F n ,L°))exp( l -J2s:(ln,L 



O exp 



a l3 (Li - Lf) (L 3 - Ll 



exp 



-^(h-L*) 



(44) 



2hK " 

Here e runs over the values {+1,-1,0} for each simplex of the triangulation, and we 
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have defined: 

(45a) 

(45b) 
(45c) 

Before applying the stationary phase approximation for small h to each of these terms, 
we need to make explicit the hidden dependence on h in the scaling^ of P°. Since h 
appears when we translate spins to areas with the LQG equivalence A = fiK^j, the 
scaling with h of P° is inverse to its scaling with j in the large spin limit: 

P a ± ~ /T m P£ (46a) 
P CT ~ ft ^Pq , < d < m , (46b) 

where P° is independent of h in the limit under consideration. 

Let us examine the numerator of (144|) . It is a sum of 3 fe terms - k being the total 
number of simplices in the triangulation - where each term corresponds to an assignment 
e = {e CT } to each simplex and has the general form 



h~ q j J\ei n AL m fS^U)0{L m ) exp 



1 

~^Qei}ni Li) 



(47) 



The power q equals the sum of the powers m or d coming from each simplex, according 
to the assignment of signs for the term under consideration, whereas f e (l n ,L m ) is given 
by pL times a product over simplices of P" and the phase function is 

Q e (l n , U) = S e (ln, U) + 7^; Oij (Li - L°) (Lj - Lfj - (U - L°) . (48) 

The "action" S e (l n , Li) has been defined as the sum over simplices of S°. For each term 
of the form (1471) . we can now apply the stationary phase expansion for small fi (which 
as we know is physically equivalent to large boundary lengths compared to lp). Just as 
in the analogous calculation in subsection I2.2[ the imaginary part of Q e is minimized 
when the boundary lengths attain the classical values peaked on by the boundary state, 
Li = L®. This selects a region D in the configuration space {l n ,Li}, and the question 
is now whether the real part of Q has an extremum in this region. The location of such 
an extremum would be given by the solution to the equations 

^-S e {ln,L°J= . (49a) 



dl n 
d 

Se{l>ni Li) 



dLi 



= K° , (49b) 

Li=L° 



15 In principle jl could also scale with h, but this scaling is irrelevant since it is an overall factor that 
cancels when dividing by Z q . 
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These equations are analogous to (I29a|l29bl) . Equation (I49al) is a Dirichlet boundary 
value problem for the action S e (l n , On the other hand, equation (149bl) is a Neumann 
boundary condition for the same action. For the particular case where e = +1 at each 
simplex, the set of equations f!49aH49bj ) has a consistent solution - the classical Regge 
solution determined by the boundary state - because the phase factor Kf is hand-picked 
to guarantee this consistency; but for a general term, the set of values {/°} that satisfy 
the first equation (solution to the equations of motion for an action with "wrong" signs 
at different simplices) is not going to satisfy the second one as well. Hence for a general 
term in the sum we fall into case (i) or (ii) of the three possibilities enumerated after 
( 1201) . and the contribution from the term is suppressed faster than any finite power of 
h. Notice that this happens regardless of the power q that comes with h in front of 
the term, therefore the dominance of the degenerate configurations in the limit of the 
vertex amplitude is now irrelevant; the semiclassical boundary state ensures that the 
contribution from the degenerate terms, as well as the sign reversed ones, vanishes from 
the semiclassical evaluation of boundary observables. 

the argument above makes a non-trivial assumption about the measure f e (l n ,Li) 
arising from the spin foam model: we have assumed that the measure vanishes smoothly 
on the boundary of the domain of integration. If this does not happen (i.e. if the 
measure gives non-zero weight to length values that saturate triangular inequalities) 
then the asymptotic expansion features a boundary term which is subdominant but not 
exponentially suppressed. As a result, restricting attention to terms with e = ±1 only 
at each simplex, a milder result holds: the "cosines" give contributions only at next-to- 
leading order. On the other hand, without further assumptions, the contributions coming 
from terms with D a factors could dominate and destroy the semiclassical behaviour. 

In summary, within our assumptions, the only non-suppressed term from the numera- 
tor is the one with an action exactly equal to the Regge action for the whole triangulation, 
with no contribution from sign-reversed or "degenerate" terms. The same thing happens 
in the denominator Z q , and thus the evaluation of observables matches that of conven- 
tional Regge calculus. The agreement is exact for the lowest order in the one-point and 
connected two-point functions, which are independent of the measure; it becomes exact 
at all higher orders as well if we posit the identification of measures 

f+(l n ,Li) :=/i(I n) i i )JIP^(i; ) Lf) = / 1 (J n ,L i ) . (50) 

This result is a strong indication that spin foam models might be able to recover the 
correct semiclassical limit, since the presence of the "Pcos S + D" factor at each simplex 
instead of the exponential of the Regge action does not affect the semiclassical evaluation 
of observables. 
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4 Discussion 



The main aim of this paper has been to make clearer and cleaner the application of 
the boundary state formalism for obtaining the semiclassical regime of nonperturbative 
quantum-gravitational theories. We have focused on Regge calculus because, as a version 
of general relativity with a finite number of degrees of freedom and a better-defined and 
relatively well-understood quantum theory it provides an ideal testing ground for the 
formalism. Moreover, Regge calculus is intimately related to spin foam models both in 
the way they are defined (from a simplicial discretization of BF theory with constraints) 
and in their conjectured semiclassical limit (through the relation of the vertex amplitude 
to the cosine of the Regge simplicial action). The main results of the paper are simply 
summarized as: 

• the semiclassical boundary state for Regge calculus successfully selects the semi- 
classical regime for the whole path integral, thereby bridging the gap between 
nonperturbative and perturbative theory for a general triangulation; 

• when at each simplex we replace the exponential of the Regge action by the cosine 
which is conjectured to come from the relevant limit of spin foam models, the results 
are unchanged because the non-Regge terms are suppressed from the semiclassical 
regime. 

The result we obtain for the connected two-point function of boundary edge lengths, a 
key object related to the simplicial graviton propagator, is 

( Li Lj ) q -{h) q { Lj ) q = h + ^i) _1 , (51) 

with S H being the Hamilton function of the configuration and a the correlation matrix 
of the boundary state. 

The approach we have followed mirrors the "graviton propagator" calculations of 
[SI [9] . The key technique in our analysis is to do the integral on internal and on bound- 
ary variables at once. The integral contains a Gaussian boundary state peaked on a 
classical configuration and we take the limit in which the relative dispersion of intrinsic 
and extrinsic geometry in the Gaussian vanish. The parameter encoding this "semiclas- 
sicality" is the ratio of the typical length scale of the classical solution peaked upon to 
the Planck scale. The stationary phase approximation, first done formally as an h — > 
limit, is revealed to be in fact an expansion in inverse powers of this large parameter. 

Having identified the parameter encoding semiclassicality allows to develop a semi- 
classical perturbation theory. The results presented in this paper are restricted to one- 
and two-point functions at the leading order, but the formalism allows us to compute 
both n-point functions and quantum corrections to the classical behaviour. The quantum 
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corrections are organized in a loop expansion akin to the one used in ordinary quantum 
field theory with a background. The difference here is that the "background" geome- 
try emerges from the nonperturbative theory only in a specific regime selected by the 
boundary state. 

Given the fact that the semiclassical boundary state plays such a key role in this 
approach, one can naturally ask whether it is possible to derive it from first principles and 
especially derive its free parameters a^. A comparison of the conventional path integral 
over all lengths (fTTil and our version that integrates only the interior and boundary of a 
region ( fT2l) suggests that the boundary quantum state can be defined as an integral over 
exterior variables: 



A possibility of this kind is discussed for instance in [51]. From this point of view, a 
choice of \l/ is equivalent to a choice of asymptotic boundary conditions in the full path 
integral. However a difficulty can be easily identified in this attempt: the a^- coming 
from the external integration is purely imaginary, therefore it cannot have a positive 
definite real part. This fact suggests that, even when considering the external integral, 
an asymptotic boundary state is needed 1 *! A standard trick to obtain the vacuum state in 
conventional quantum mechanics and field theory is to perform a Wick rotation. However 
the implementation of a meaningful analogue of this procedure (and even the existence 
of one) in a general-boundary setting is still unclean 17 !. 

A more promising route might be to select the boundary state representing the semi- 
classical regime in an a posteriori way, by requiring that the results match those of 
conventional perturbative Regge theory (in which only small fluctuations around a given 
classical solution for the whole manifold are quantized and summed over; see [60|, [6T ] [62] ) . 
Clearly this vindicates the general Gaussian expression to the lowest order, and the free 
parameters in this parametrization (one for each pair of boundary edges Li, Lj) can 
be established a posteriori by matching the correlations between each pair of bound- 
ary edges with the conventional perturbative expression for it. A nontrivial consistency 
check of this matching can be found if we consider deformations of the boundary. Given 

16 As a side remark we would like to recall that, in the standard approach to perturbative quantum 
field theory where states are associated to spatial sections, the vacuum state a t a given time t is also 
a Gaussian functional with a real positive definite a. It can be obtained from the partition function 
Zq by integrating out from time t to — oo; however, unless a Wick rotation is implemented, this would 
give an imaginary a. This shows that even in ordinary quantum field theory the partition function 
Zq contains a boundary vacuum state with real a (defined asymptotically at t — ±00). This issue is 
discussed in chapter 9 of Weinberg's [70] and is the origin of the +ie prescription. 

17 In [T31 U31 HS1 USl HZ] the boundary vacuum state is found starting with a Gaussian ansatz and then 
imposing a set of consistency conditions. 
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a Gaussian boundary state for a given boundary including edges I and J, if we choose a 
different boundary which also includes these two edges we can find (by integrating out 
from (1521) the edges external to one boundary and not to the other one) how the new co- 
efficient a'jj in the new boundary state is related to the old one ajj. On the other hand, 
the requirement that the correlation between two edges is invariant under the change of 
boundary (which must be true to match conventional perturbative theory, which knows 
nothing about boundaries) gives a separate constraint on the variation of the a matrix, 
and the consistency of the two constraints is a nontrivial check on the validity of the 
approach and its ability to recover conventional perturbative theory correctly. We hope 
to discuss these matters more fully in a future paper. 

Moving on to the results of Section [2J2], we believe that they may help clarify an impor- 
tant aspect of the "gravitons from spinfoams" calculation: namely Rovelli's suggestion 
that few spinfoam vertices can provide a good approximation for low energy processes 
[7T] . We find that, in the flat space case for quantum Regge calculus, connected two- 
point correlation functions are, at leading order, independent of the triangulation. This 
means that, given a triangulation of the boundary, a minimal triangulation of the bulk 
(one with no internal vertices, so that only a number of bulk diagonals is present) is 
enough to obtain the exact result. Triangulation independence is expected to be lost at 
next-to-leading order; as a result we would have that quantum corrections to simplicial 
flat space correspond to unfreezing the degrees of freedom associated to the internal 
vertices of the triangulation. 

Quantum Regge calculus provides an ideal arena for testing the ideas introduced in 
[H [9]. Moreover, the results presented in this paper (and in particular in section E]) 
may turn out to be directly relevant for the analysis of the semiclassical regime of the 
new spin foam models [3U1 ED E21 ESI Ell beyond the single vertex level, if a number 
of gaps are filled. The first one is obviously whether their large spin limit is in fact 
related to the cosine Regge expression for each simplex as postulated in expression (|42|) . 
There are strong indications that this is likely to be true. On one side, recent numerical 
and analytical work on the so-called "flipped" vertex [721 [73] indicates that it has a well- 
behaved semiclassical regime with a geometrical interpretation; on the other, a stationary 
phase analysis of the Freidel-Krasnov (FK) model [33] after writing it in path integral 
form [711 EH] finds not only that the semiclassical limit of the vertex is the cosine of 
the Regge action, but also that the mismatch between area and length variables is side- 
stepped since area configurations that do not correspond to simplicial geometries drop 
out from the semiclassical limit. This last development is especially encouraging since the 
variable mismatch [76] and the presence of discontinuous boundary geometries [23l 164"] 
is another general problem the spin foam formalism seems to face when trying to obtain 
a semiclassical limit with a clear geometrical interpretation. 
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The recent formulation of the FK model as a path integral [74] also opens the door 
for the possibility of doing the stationary phase analysis with a boundary state in a 
purely spin foam context. Hopefully this would vindicate the implicit assumption in our 
"hybrid" expression (j4"2"l) that the semiclassical "cos S + D" expression for the vertex 
amplitude can be taken at each internal vertex of the two-complex, even if the state sum 
is nonperturbative and the semiclassical regime is enforced with a boundary state that 
peaks explicitly only the boundary variables. It seems likely that this could be done 
by adding a suitable boundary state to the formal stationary phase analysis of the FK 
model mentioned in the previous paragraph. 

We believe that the results of section [3] about the dropping out of the all the non- 
Regge terms coming from the conjectured limit of the spin foam vertex amplitude is very 
encouraging news for the spin foam research program. The "cosine problem" has been 
a potential obstacle for the semiclassical limit of spin foams from the very beginning 
of this research program, as can be seen already in Ponzano and Regge's 1968 paper 
[55] , where they say referring to the "non-Regge" terms arising from the large spin limit 
in 3d gravity: 11 The other terms, other than the positive frequency part, are related to 
different orientations of the tetrahedra Tj and have a similar interpretation, although 
their precise meaning is still unclear. It is plausible that in the transition to a smooth 
manifold M. they will give no essential contribution to the final result" . Our results show 
that if the large spin limit of a spin foam state sum is a Regge-like path integral with the 
cosine factor at each simplex, Ponzano and Regge's hope is realized and the non-positive 
frequency terms give no contribution to the final result, once the desired semiclassical 
regime is selected with an appropriate boundary state. 
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